m 



Collisions, magnetization, and transport coefficients 
in the lower solar atmosphere 

J. Vranjes^ and P. S. Krstic^ 

^Institute of Physics Belgrade, Pregrevica 118, 11080 Zeniun, Serbia 
^Joint Institute of Computational Sciences, University of Tennessee, Oak Ridge, TN 

37831-6173, USA 

Abstract: The lower solar atmosphere is an intrinsically multi-component and colli- 

^ sional environment with electron and proton collision frequencies in the range 10^ — 10^° 

5-1 Hz, which may be considerably higher than the gyro-frequencies for both species. Colli- 

<^ sions between different species are altitude dependent because of the variation in density 

lO and temperature of all species. 

We aim to provide a reliable quantitative set of data for collision frequencies, magne- 
f^ tization, viscosity, and thermal conductivity for the most important species in the lower 

C/^ solar atmosphere. Having such data at hand is essential for any modeling that is aimed 

at describing realistic properties of the considered environment. 

The relevant elastic and charge transfer cross sections in the considered range of colli- 
sion energies are now accepted by the scientific community as known with unprecedented 
^, accuracy for the most important species that may be found in the lower solar atmo- 

sphere. These were previously calculated using a quantum-mechanical approach and were 
t>" validated by laboratory measurements. Only with reliable collision data one can obtain 

T-H accurate values for collision frequencies and coefficients of viscosity and thermal conduc- 

^ tivity. 

■^ We describe the altitude dependence of the parameters and the different physics of 

O 

(y-) collisions between charged species, and between charged and neutrals species. Regions 

of dominance of each type of collisions are clearly identified. We determine the layers 
> 
•'^ within which either electrons or ions or both are unmagnetized. Protons are shown to 

;_( be un-magnetized in the lower atmosphere in a layer that is at least 1000 km thick even 

for a kilo-Gauss magnetic field that decreases exponentially with altitude. In these layers 

the dynamics of charged species cannot be affected by the magnetic field, and this fact is 

used in our modeling. Viscosity and thermal conductivity coefficients are calculated for 

layers where ions are unmagnetized. We compare viscosity and friction and determine 

the regions of dominance of each of the phenomena. 

We provide the most reliable quantitative values for most important parameters in 

the lower solar atmosphere to be used in analytical modeling and numerical simulations 

of various phenomena such as waves, transport and magnetization of particles, and the 

triggering mechanism of coronal mass ejections. 
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1 Introduction 

There has recently been a considerable shift of focus of solar researchers from ideal and 
collision- less toward coUisional phenomena in the solar atmosphere (e.g., Arber et al. [2], 
Pandey and Wardle |29], Soler et al. [32], Barcelo et al. |3], Zaqarashvili et al. [S], |i5] . 
Khomenko and Collados [2T]). This is understandable for lower atmosphere layers because 
this is a partially ionized medium with several species whose dynamics is heavily coupled 
and affected by mutual collisions and by collisions between similar particles (Vranjes and 
Poedts [3H], Vranjes et al. [ID], H^). In some layers of the photosphere the ion- neutral 
and electron-neutral collision frequencies are approximately 10^ and 10^° Hz, respectively. 
This makes both ions and electrons in these layers very weakly magnetized or completely 
un-magnetized. 

In recent multi-component models in the literature, the medium has typically been 
treated as if it consisted of two components, neutrals and 'plasma'. These models con- 
sequently neglected differences between electron and ion dynamics, similarly to ordinary 
magneto-hydrodynamics theory. However, in the lower solar atmosphere it is very difficult 
to justify this approach, as we show below. One reason for this are the different collision 
frequencies of electrons and ions, which among other effects imply different magnetization 
of these two species, consequently producing quite different effects of the magnetic field 
on particle dynamics. 

The key components in the lower solar atmosphere are identified in the present work. 
We also show their collision cross sections and collision frequencies as a function of al- 
titude. Using these results, the altitude-dependent magnetization of electrons and ions, 
the viscosity tensor, and the thermal conductivity vector components are calculated. The 
altitude-dependence of the parameters is presented graphically. 

2 Key ingredients 

The lower solar atmosphere is an essentially multi-component system that includes elec- 
trons, protons, neutral hydrogen H and helium He atoms, helium ions He"*" or/and He^"*" 
etc. The dissociation energy of molecular hydrogen H2 is ~ 4.5 eV and its presence may 
be expected throughout the lower solar atmosphere. However, recent observations (Jaeg- 
gli et al. [IB]) showed that the regions with the highest amount of molecular hydrogen 
are sunspots where the total fraction of H2 molecules may reach a few percent only. On 
the other hand, the ionization energy of the hydrogen molecule, 15.603 eV, exceeds the 
ionization energy of atomic hydrogen, so the presence of the ionized hydrogen molecule 
is most likely negligible. Having in mind all possible collisional combinations of various 
particles, it is essential to identify the most important of their interactions before making 



any attempt to obtain realistic models for wave damping, diffusion, transport, magneti- 
zation of particles, etc. The collision cross sections can vary with the energy of colliding 
particles, which in the solar atmosphere is equivalent to the variation of the temperature 
with altitude. 

Detailed studies available in the literature (Krstic and Schultz [23], Glassgold et al. 
[16], Schultz and Krstic [31], Jamieson et al. [19] etc.) show large differences between 
the cross sections describing coUisional phenomena such as elastic scattering, momentum 
transfer, viscosity, spin exchange, and charge exchange. Among these the most prevalent is 
typically the elastic scattering. Physically, this cross section should be taken into account 
in estimating the magnetization. On the other hand, the cross section for momentum 
transfer should be used in calculating friction; this type of cross section turns out to 
be usually lower than the elastic scattering. Knowing these details may be essential in 
realistic modeling of the lower solar atmosphere. 

One also needs to include some inelastic collisions, which clearly may play a decisive 
role in the partially ionized solar atmosphere (Vranjes and Poedts [SS]), which makes 
the whole analysis significantly more complicated. These processes are the radiative 
recombination (of the type A'^ + e — )■ A + hu), three-body recombination (of the type A+ + 
e+e — )• A+e), dissociative recombination (e.g. of the type Ag'+e — )■ A+A*, where A* is an 
excited atom), ionization by electron impact (the process of the type A + e ^ A^ + e + e), 
the ion-atom (or molecule) interchange reactions of the type A'^ + BC — )■ (AB)^ + C , etc. 
The presence of inelastic collisions implies modifications of equations through appropriate 
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Figure 1: Integral cross sections apn for proton {p) collisions with neutral hydrogen atoms 
H for quantum-mechanically indistinguishable nuclei of the projectile and target particles, 
following Krstic and Schultz [23] ■ Here and throughout the text 1 a.u. = 2.8 
le1/~ 11604 K. 
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source/sink terms that appear in the continuity equation, in the momentum and in the 
energy equations. Only collisional phenomena from the first group mentioned above, i.e.. 



elastic and charge transfer processes, are discussed here. Even then, as will become evident 
in the following sections, we still encounter a plethora of processes that could further be 
classified by priority and dominance to perform any meaningful analysis. There exists a 
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Figure 2: Integral cross sections apne for proton collisions with neutral helium He. 

numerous literature dealing with collisional plasmas in the laboratory environment. The 
laboratory knowledge is in principle directly applicable to the relevant phenomena in the 
solar atmosphere. However, we find that i) the theoretical and experimental results are 
scattered to a large extent, and, ii) for the purpose of solar plasma studies, the cross 
sections for a specific range of parameters of the solar plasma are missing. For this 
reason we present cross sections for possible key-players in the lower solar atmosphere in 
the following subsection, with the energies of colliding particles in the range appropriate 
for the plasma below the transition region. These results, combined with the existing 
models of the altitude-dependent densities and temperature of the plasma species, are the 
essential and sufficient ingredients for calculating collision frequencies that are necessary 
for modeling. 



2.1 Cross sections for charge- neutral and neutral-neutral colli- 
sions 



Figs. [T]j4] described in this section provide cross sections for all most important ingredi- 
ent species in the lower atmosphere, given in terms of the energy (temperature) of the 
colliding species, which makes them directly applicable to the varying temperature with 
the altitude. The cross section profiles in Figs. [l]J3] are based on works of Krstic and 
Schultz |23], [21], [25]. These authors derived both differential and integral cross sections 
for elastic scattering. The cross sections for the momentum transfer and viscosity are the 
first and second moments of these [see also Dalgarno et al. [TT], Hasted [TTj, Brandsen 
[S], Makabe and Petrovic [2S], Schunk and Nagy [33]] • In the first (momentum transfer) 



the differential elastic cross section is weighted by a scattering angle 6 dependent term, 
1 — cos 6*, while in the second (viscosity) it is weighted by sin^ 6. These weighting factors 
describe different features associated with momentum transfer and viscosity. For the vis- 
cosity this emphasizes the scattering at an angle 7r/2 and de-emphasizes the forward and 
backward ones. Scattering to such large angles is very effective in equalizing energies of 
the colliding particles. This is seen from the expressions for energies of the two particles 
before {Ei, E2) and after collisions {E^, E2) in the laboratory frame, expressed through 
their total energy E in the laboratory frame and scattering angle 9 in the CM frame: 
E^ = E{1 + cos6')/2, E2 = E{1 — cos0)/2. Hence, such large-angle scattering tends to 
reduce both viscosity and conductivity. On the other hand, the factor 1 — cos 6' in the 
momentum transfer cross section emphasizes the backward-scattering angles, and this 
cross section determines the average momentum lost in collisions. Note also that charge 
exchange is a backward-scattering process. Many more details on these cross sections can 
be found in Krstic and Schultz [23], [25], currently the most accurate cross sections for 
elastic processes and resonant charge transfer. The energy range in Figs.[l]j3]in the center 
of mass (CM) of colliding particles is 0.1 — 5 eV (bottom x-axis), while in the laboratory 
(plasma) frame the energy range is given by the top x-axis using the transformation for- 
mula Eiab = EcMij^i + f^2)/i^2, where m2 is the mass of the target particle, and 1 a.u. 
= 2.8 ■ 10-21 m^, le\/ ~ 11604 K. 

We start with the cross sections for proton collisions with neutral hydrogen (p-H), 
shown in Fig. [ij They are based on quantum-mechanical indistinguishability of the pro- 
jectile and target nuclei (Krstic and Schultz [23] )■ The elastic scattering curve from Fig. [I] 
(line 1) was used to calculate magnetization, i.e., for the ratio of the collision frequency 
and the gyro-frequency of protons. It is the sum of the pure elastic scattering cross sec- 
tion and the charge transfer cross section. In estimating the friction caused by neutral 
hydrogen, it is appropriate to use the momentum transfer cross section (line 2). The 
amplitude oscillations of the cross sections are the consequence of quantum effects, which 
are present only at lowest collision energies (in the present study this means throughout 
the photosphere and chromosphere). 

We stress that when one approximates the distinguishable particles, the elastic scat- 
tering cross sections and their higher momenta (momentum transfer and viscosity) are 
lower. This is because in this case we assume that we can distinguish between elastically 
and charge-transfer scattered protons, resulting in a separate treatment of these processes. 
When these particles cannot be distinguished, as is the case at lowest energies (lower than 
1 eV), the elastic cross sections of the indistinguishable particles and their moments are 
the coherent sum of the processes, elastic and charge transfer. For indistinguishable nu- 
clei, the integral elastic cross section together with the charge transfer at 0.2, 0.5, 1 eV 
in the CM frame is, 788.660,679.534, and 582.292 a.u., while the distinguishable nuclei 



model yields 598.068, 506.701, and 419.739 a.u. Similar holds for the momentum transfer, 
while the charge transfer cross section is practically the same for both models. One has to 
have in mind these differences and the differences in the physical definitions of these cross 
sections, to avoid twice counting elastic and charge transfer cross sections: the elastic 
cross section of the indistinguishable particle and their moments in Fig. [T] already contain 
both processes. 
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Figure 3: Integral cross section ann for collisions between neutral hydrogen atoms H for 
quantum-mechanically indistinguishable nuclei. 



In Fig. |2] the three lines describe the collisions between protons and neutral helium 
atoms. Observe that the momentum transfer and viscosity lines are below the line for 
elastic scattering by a factor 4-5. This is because the momentum transfer presents the 
differential cross sections in the backward-scattering directions, while the dominant con- 
tribution in the elastic cross sections comes form the forward-scattering angles, which 
dominate the differential elastic cross sections. Therefore, for the proton dynamics the 
presence of neutral helium may be more important for estimating magnetization than for 
the momentum loss caused by friction or by viscosity (as compared to those that come 
from proton self-collisions or interaction with hydrogen, see more in Sec. |6]). 

Fig. [3]gives collision cross sections between neutral hydrogen atoms according to Krstic 
and Shultz [23]. The momentum transfer and elastic collisions curves coincide, and line 
1 in the two figures describes the most dominant interaction for the solar atmosphere. 
Note that it includes both direct and recoil scattering as a direct consequence of the 
indistinguishability of particles, and the same holds for the viscosity cross section. Due 
to these reasons the presented values are twice as high as the classic values obtained from 
the model of distinguishable particles. 

The new type of cross section that appears in Fig. [3] is the spin exchange cross section, 
which describes collisions in which electrons (from the two colliding atoms) with different 
spin orientation are exchanged. Processes of this type are the reason for the cooling 




2 4 6 

electron energy [eV] 



10 



Figure 4: Cross section for elastic scattering of electrons on neutral hydrogen atoms H 
and neutral helium atoms He, in terms of the electron energy. 



phenomena in the upper atmosphere, interplanetary space, and galactic HI regions (Bates 
m, Purcell and Field ^30], Dalgarno [12]). 

Observe the difference between the lines 1 in Fig. [T] and in Fig. [3j One obvious reason 
for this is the charge transfer cross section, which is contained in line 1 in Fig. [ij This can 
be subtracted to obtain the pure elastic scattering value <7ei,pH- The cross section obtained 
in this way, at high energies, tends smoothly toward the corresponding elastic cross section 
obtained from the classic model of distinguishable particles. At the low energies of interest 
here, 0.5 and 1 eV (in CM frame), this yields aeLpn = 507.333,420.038 a.u. as the pure 
elastic scattering cross section for p-H collisions. For the H-H collisions in Fig. |3] the 
corresponding values are lower by about 137 and 81 a.u. for the two energies, respectively. 
The difference that still remains (i.e., greater p-H than H-H cross section) should be 
attributed to the fact that proton collision with H atoms causes charge polarization on 
the neutral atom [see Chen and Chang [TUj and in Vranjes et al. 137]] • This means that 
to some extent p-H collisions involve features of the Coulomb interaction. The physics 
behind this is as follows: a point charge go placed at some distance from an atom with 
the radius a that has a point positive charge q in the core and uniform negative cloud —q 
around it, will cause displacement of the initially uniform cloud charge. If the external 
point charge go is an electron, this displacement is in the direction away from the electron 
position. If the external point charge go is a positive ion, this displacement of the cloud 
charge will be toward the ion and the force will again be attractive. The energy from this 
attractive Coulomb interaction according to Chen and Chang [10] is 
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Hence, collisions of charged species with neutrals indeed involve a Coulomb-type inter- 
action, which affects collisions at very short distances, see also Dalgarno et al. [11] and 



McDowell and Coleman 

Finally, we give plot in Fig. |4] with three lines for the cross sections for the electron 
scattering on the two most important atoms hydrogen H and neutral helium atoms He in 
the electron energy range 0.1 — 10 eV. The lines represent some mean values from many 
references, see for example Bedersen and Kieffer [5] and references cited therein. The 
possible uncertainty is almost 10-25 percent at low energies. This also agrees with some 
other sources, for instance Erode [7], Brackmann et al. [6], Kieffer [22], Mitchner and 
Kruger [2H], Tawara et al. [35], and Fortov et al. [T5]. Hence, although some uncertainty 
for electron cross sections exists, it is not substantial. Clearly, at the low energies of 
interest for photosphere and chromosphere the most probable are electron collisions with 
atomic hydrogen H. We investigate the electron Coulomb collisions in the following section. 

3 Electron collision frequencies 

To describe the collisions between charged particles we use the following expression for 
the collision frequency between the charged species b and a following Spitzer [3lj and 
Vranjes et al. |13]: 
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As is well known, the Coulomb logarithm Lfya (introduced by Spitzer) describes the 

cumulative effect of numerous small angle deflections that are intrinsic to Coulomb-type 

collisions. 

These expressions are used to calculate collisions for electrons for the parameters 

(density and temperature) that vary with altitude. The results are presented in Fig. |5} 

To incorporate the variation of the parameters, here and throughout the text we use the 

values given in Table C in Fontenla et al. |14j . 

For electron collisions with neutrals we read the cross sections cr(x) from Fig. Ill and 

then calculate the collision frequency using the expression 



^ej (^ej 



[x)nj{x)vTe{x). (3) 



The small differences of the parameters (density and temperature) given in the reference 
above, as compared with some other models of the lower solar atmosphere from the same 
authors or others, are of no particular importance for the general picture that is obtained. 
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Figure 5: Electron collisions with the altitude. For comparison, the thin line gives the 
electron gyro-frequency Qe for the starting value of the magnetic field Bq = 0.1 T. 



The same holds for the expression h>ej used here in comparison with some modifications 
of it that may be seen in the literature. For example, the thermal speed we use is 
without any numerical parameter, as for the mean velocity v = [8kT / (7imab)Y^'^ , "^ab = 
TTiafnb/i'tTT'a + ^^ifc), which is sometimes used in the literature. It is easily seen that in the 
most drastic case, for example when a = b, this increases our thermal speed by a factor 2.2 
only. Similar numerical parameters appear in the most probable speed v = {2kT /mY^'^^ 
and in the root-mean-square velocity (-O^)^/^ = (SKT/mY'''^. However, these modifications 
are not substantial in view of our much more accurate cross sections as compared with 
those that are typically used in the solar plasma literature (see comments in Sec. ffl. 
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Figure 6: Coulomb collisions of protons in terms of altitude. 

The electron collision frequency was checked also for e-He"*" collisions. For the given 
altitude range the maximum e-He"*" collision frequency is at about x = 2000 km altitude, 
but it is only 0.02 [in the same units as in Fig. [s] and is therefore completely negligible. 
After checking for the electron collision frequency with He"*"^ ions we found out that it 



was even lower, at least by one order of magnitude. 

The electron collisions with neutral helium He are even lower than the dominant 
collisions in Fig. [sj For example, at x = we have Uene = 2.2 ■ 10^ Hz, which is almost 
two orders below u^h, and it remains well below in the whole region. 

Fig. |5] suggests that in the region — 850 km the electrons' collisions with atomic 
hydrogen are by at least two orders of magnitude more frequent than the electron-electron 
collisions. Above 850 km e-p collisions (and the associated friction) should be more 
important than both e-e collisions and electron collisions with neutrals. Below the level 
denoted by x = the Coulomb collisions become more dominant. 
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Figure 7: Collision frequency for protons colliding with the neutral atoms H and He. 



According to Table C in Fontenla et al. [14j, at 850 km altitude the proton and H 
number densities are about 10^^ m^^ and 10^° m~^, respectively, i.e., neutral hydrogen 
atoms are about 1000 times more abundant, yet e-p collisions are already more frequent. 
A similar situation is observed in the interval between x = —100 km and a; = 0. Here and 
throughout the negative altitude denotes the value below the referent level x = 0, i.e., 
toward the center of the Sun. This all confirms the well-known fact that the Coulomb 
collisions have a much larger cross section and are more frequent even in rather weakly 
ionized plasmas (see in Ratcliffe [IS], Vranjes and Poedts [IJ)- These facts are frequently 
overlooked in the literature. 



4 Proton collision frequencies 



Proton collision frequencies Upj where j includes protons as well as other relevant charged 
or neutral species are presented in Figs. |6] and [7| Here, for proton-neutral collisions 
we have Upj = apjUjVTi and apj is given in Figs, [l] and [2] in lines 1 and 2. In Fig. [6] the 
Coulomb p-p and p-He+ collision frequencies are presented using Eq. ^ . Although p-He"*" 
collisions are obviously much less frequent, their actual importance may be understood 
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only by comparing the friction and viscosity terms in the momentum equation, see Sec. 

El 

In Fig. [tI the cross section a{x) = a(T(x)) is obtained from Figs. [I] and [2] for the 
energies in the plasma frame (i.e., those given by the top x-axes), and the density and 
temperature (energy) from Table C in Fontenla et al. [11] . 

Comparing Figs. |6] and [7} clearly proton collisions with neutral hydrogen are by far 
the most dominant in the given range (which starts from x = —100 km) up to the altitude 
of about 1350 km (there is a difference of almost four orders of magnitude for collisions in 
certain lower layers). Above 1350 km the amount of neutrals is sufficiently reduced so that 
p-p collisions become dominant. The highest proton collision frequencies at x = —100 km 
read Up^ = 2 ■ 10^ Hz, Upp = 7 ■ 10^ Hz, Up^e = 8-10^ Hz. 

Above 1900 km proton collisions with He"*" are more frequent than p-H, which is even 
more true for p-He collisions. For example, at a; = 2017 km we have i'pHe+ = 3.3 ■ 10^ 
Hz [see Fig. |6], while Up^ = 1.2 ■ 10^ Hz and z/p^e = 62 Hz [see Fig. [7J, and higher up 
this difference increases. Accordingly, above this layer proton friction with (any) neutral 
atoms is negligible. 

5 Electron and proton magnetization 

To estimate the magnetization, the thin line in Fig. [5] gives the electron gyro-frequency 
fie(a^) = ^Bq{x) / rUf^ for a starting value of the magnetic field Bq = Q.IT which approx- 
imately changes exponentially with the altitude as exp[— a;/(2/i)], h = 125 km, following 
the thin-flux tube model and the usual pressure balance conditions. In the region below 
X = the Coulomb collisions (e-p and e-e) become more dominant [see Fig. Is], and clearly 
the collision frequency in that area is higher than the electron gyro-frequency even for 
Bq = 0.1 T. Hence, in this layer Ve/^e > 1 and electron dynamics should not be influenced 
by the magnetic field (see also Vranjes et al. [lO]). 

The layer without magnetization is much wider for protons and other ions. The 
corresponding line for protons is given in Figs. [6] and [7| For the same exceptionally strong 
starting field i?o = 0.1 T we have fli ~ 9.6 MHz, which changes with altitude in such 
a way that it remains below Up^ up to at least 1000-1200 km. Above this altitude the 
magnetic canopy is formed and the field changes less rapidly; it is therefore expected 
that above this altitude the proffie for rii{x) is less steep. From Figs. |6]and[7]it is seen 
that ions remain un-magnetized within a layer of unknown width for x < 0. Accordingly, 
because no much stronger magnetic field can be expected elsewhere in photosphere, we 
can conclude that there exists a layer throughout the photosphere which is at least 1000 
km thick (most likely it is even thicker) within which protons remain un-magnetized in 
absolute sense. 
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Geometry and magnitude of the magnetic field vary both horizontally and vertically. 
Therefore the width of the layer within which protons are unmagnetized might be expected 
to be much wider in regions with a considerably weaker field. However, assuming that the 
magnetic canopy forms at about the altitude of 1000 km (c.f. Khomenko et al. [2U]), this 
implies that protons are unmagnetized below the canopy; this holds throughout the lower 
atmosphere. In any case, with the accurate values for collision frequencies presented in 
Figs. |5|J7| the actual magnetization and the width of this layer can easily be checked for 
any given value of the magnetic field. 

We can claim with certainty that there exists a well-defined (but highly irregular 
regarding its width) altitude range within which both electrons and protons are totally 
unmagnetized; this fact should not be ignored in modeling. The mean free path of a 
particle j (the distance it covers between two consecutive colhsions, given by Xjj = Vj-j/uj) 
in these regions is far shorter than the ion gyro-radius. Hence, the dynamics of ions is 
not affected by the magnetic field in most of the photosphere and chromosphere. In 
some layers this holds for electrons too, and such an environment can support only waves 
appropriate for an ordinary gas (e.g. gravity and/or acoustic oscillations) or heavily 
damped ion-acoustic waves. 

In addition to this, according to numbers presented above, there exists an altitude 
region within which electrons are magnetized while protons (ions) are not. This makes it 
very difficult to justify so called two-component models that are found in recent studies 
which assume the medium to be composed of neutrals from one side and 'plasma' from the 
other. The term 'plasma' here refers to electrons and ions as a single fluid. Because there 
are magnetized electrons and unmagnetized ions in the regions that we clearly identified, 
we know that the dynamics of the two species perpendicular to the magnetic field or at 
large angles with respect to it becomes totally different, which precludes describing them 
with a common set of single-fiuid equations. 

6 Viscosity and thermal conductivity in unmagneti- 
zed plasma 

Because of the altitude-dependent parameters, the contributions of different species to 
viscosity and conductivity coefficients will vary in space, and the spatially dependent 
contribution of each component should be checked separately. The viscosity tensor com- 
ponents i,j for species a in a strongly collisional plasma-gas mixture with un-magnetized 
charged species are given by 

n _ Pa ( dVa,i dVgj 2 dVgA 

"'*^ Eft ^ab V 9^3 d^i 3 '^ dVk ) 
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Here, rjj^fc stands for the coordinates x, ?/, z, while f j^-^a; denotes the speed components 
along these coordinates, and summation (with the general index h) includes all species 
including the specie a itself [i.e., collisions between alike particle as well, when the terms 
in the second row in dil) clearly vanish]. From Eq. Q it can easily be seen that this is a 
symmetric tensor IIjj = IIjj, and its trace is zero, 11 jj = (with assumed summation over 
the repeating index). These are well-known features of the viscosity tensor. 

The first row in Eq. ^ describes the self-induced viscosity of the species a. The 
second row on the other hand is due to relative motion of the species a with respect to 
other species (this implies collisions between dissimilar particles). This part may play a 
key role in the initial stage of some accidental electromagnetic or electrostatic perturba- 
tions in which other (uncharged) species are at rest; this holds for the friction force and 
friction damping as well. Because of collisions between dissimilar particles, Eq. ^ differs 
considerably from the usual Navier-Stokes formula for single-component gasses. 

The components of the conductivity vector are given by 
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Similar to the viscosity, here the term in the first row in Eq. ([s]) is also due to self- 
conductivity and the remaining terms include interactions between dissimilar species. 

Eqs. (Il]), ([5]) are obtained from kinetic equation with the Bhatnagar-Gross-Krook 
(BGK) collisional integral. The Grad method is used together with the fact that the 
temperature variation in the photosphere-chromosphere layer studied here is very weak, 
it changes for about 0.3 eV only. The general transport coefficients contained in Eqs. ^, 
(p| are given in Sec. 6.1 



In more explicit form, the components of the viscosity tensor for the un-magnetized 
species a are 
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In the expressions above it is appropriate to take pg = ngKT. Hence, all species have the 
same temperature and there is no anisotropy. Both assumptions are well-justified in such 
a strongly collisional and un-magnetized lower solar atmosphere. 

6.1 Proton dynamics 

Using all previous graphs, we can now calculate the viscosity coefficients for the particular 
photospheric plasma. For protons, after identifying the leading contributors to their 
collisions in Sec. [4], we need the following coefficient for the viscosity in the first row of 

Eq. Q: 

np(x)/tT(x) 
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[in kg/(sm)]. 
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In the second row of Eq. Q we have the following four coefficients [all in units kg/m^] : 
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Figure 8: Dynamic viscosity coefficients ([6]), ([8 10) without the magnetic field for protons 
in the lower solar atmosphere. Here rjpp is plotted in units 10~^ kg/(s m) and all other 
coefficients in 10"^'' kg/m^. 



Note that to calculate fipnix) and fipnei^) "we have to use the viscosity lines (line 4 from 
Fig. IT] and line 3 from Fig. |2| respectively) and the top x-axis for the energy. 

The numerical values for the coefficients ^ and (|8 10 ) are presented in Fig. [s] The 
coefficient fXpp is not presented because from the second row in Eq. dll) it is seen that 
for a = 6 its contribution vanishes. Therefore the most dominant fiab terms should be 
checked only for the case a ^ b. From Fig. [8] it is clear that in most of the space the 
coefficient fip^ should be taken into account. However, above an altitude of about 1700 
km the viscosity between protons and helium ions becomes the most dominant, as is seen 
from the dotted line, which gives the values of fipHe+, though this may change if protons 
are magnetized; then the gyro- viscosity should be taken into account. This problem will 
be discussed elsewhere. Furthermore the viscosity that involves neutral helium fipj^e is 
clearly negligible everywhere. 

It is meaningless to directly compare the leading fiab term with r]pp because they are 
in different units. Instead, one must compare the complete corresponding viscosity terms 
from the first and second row in Eq. Q, which is approximately 
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For waves with the wave number k we have L^ ~ 1/k. Because the ratio rjpp/fipH changes 
with altitude for many orders of the magnitude, the relative contribution of the two terms 
to the viscosity will drastically change. For example, at x = we have rjpp/fipH = 0.22 
[m^/s] and at x = 1580 km rjpp/fipH = 3.2 ■ 10"^ [m^/s]. The speed difference 6vi may be 
time dependent, for instance for waves that first affect some charged species a while the 
un-charged species b is set into motion only after some coUisional time, which will also 



affect the ratio (11) 
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Without relative motion between protons and other species, the proton viscosity is 
strictly self-induced. With relative motion the situation becomes much more complicated 
and it is not so obvious which terms are more dominant. This is true in particular in view 



of ratio (16), which involves friction (see below). Hence, to be on a safe ground, one should 
keep the two leading viscosity terms discussed above together with the corresponding 
friction force terms. 

A similar analysis is now used for the conductivity vector (pi). For unperturbed sit- 
uations the second row in Eq. ^ can usually be neglected because i) collisions couple 
separate species and they are expected to move together, and ii) in the presence of fre- 
quent collisions thermalization is very effective, so the temperatures are equal. This does 
not necessarily hold for some accidental electrostatic or electromagnetic perturbations. 
For example, electrostatic ion-acoustic-type perturbations in the given environment must 
involve ion temperature perturbations (because the wave phase speed is on the same order 
as the ion thermal speed), and this happens on the background of initially static neutrals. 
Therefore the second row in ([s]) should be kept. For protons, keeping the most relevant 
terms as above, we have the x-component of the conductivity vector 
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In view of Eq. (Is]), the thermal conductivity coefficients which appear here are 
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Here, similar to Eqs. ([6 10), all parameters are altitude dependent, which affects the ther- 



mal conductivity coefficients (12 15), whose altitude dependence is presented in Fig. |9jin 
units 10^° (sm)~^ for Upp, and 10"^" kg/m^ for Xpb coefficients. Similar to viscosity, regard- 
ing Xpb coefficients, here again proton interaction with hydrogen atoms is most dominant 
below 1700 km and only Xph should be kept in the part that describes the interaction 
between different species in Eq. ([5|. Above this altitude one should keep XpHe+ only. 
These conclusions hold as long as the ions are un-magnetized (gyro- viscosity excluded). 
Identifying these most dominant terms can significantly simplify the derivations. 
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Figure 9: Proton thermal conductivity coefficients (12 15) with altitude; Kpp is plotted in 
(sm)~^, and Xpj iii kg/m^- 

6.2 Friction vs. viscosity. Vanishing friction 

Comparing the contribution to the momentum equation of the viscosity term from the 
second row in Eq. ^ and the usual friction force yields the following approximate dimen- 
sionless ratio of these terms for the species a: 
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Ff 0-m,ab Hb ^a 



(16) 



Here, we estimate the interaction between the species a and another single species 6, and 
o'v,abl(^m,ab < 1 is the ratio of the cross sections for viscosity and momentum transfer [which 
are typically different, see Figs, ^p], and k^^ = A is the characteristic scale length for the 
speed difference (e.g., wave length for wave analysis), which appears from V ■ ^a,ij in the 
momentum equation. It is clear that Ri^a can have any value [e.g., because the difference 
Va — Vb may be time dependent for electromagnetic or electrostatic wave phenomena that 
affect the charged species a ffist]. Therefore it cannot be justified to neglect the second 
row in Eq. ^ in the wave analysis. Assuming that a represents charged species and b 
some neutral one, for electromagnetic or electrostatic perturbations species b is initially 
at rest so that Ri^a < 1 if 

X> '—. 

l^b ^ab <^m,ab 

For these wavelengths the contribution of the second row in the viscosity in Eq. ^ is 
negligible, but this holds in principle for the initial regime only. 

Comparing now the viscosity contribution to the momentum equation due to the ffist 
row in Eq. (|4]), with the friction force between a and b for the initial stage with v^ = 0, 
yields 
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Table 1: Common velocity and velocity relaxation time tc for protons (with arbitrary 
starting speed Vpo) and hydrogen (with Vno = 0) for several altitudes. 



X [km 


650 1180 2043 


Vc 
tc [s 


5.1 ■ 10-%o 2.2 ■ IQ-^Vpo 9.5 ■ lO'^Vpo 0.37t;po 
5.6-10-9 1.2-10-6 7.6-10-5 5-10-3 



The friction force dominates if 



?;2 



The index m here denotes the collision frequency calculated with the appropriate cross 
section for momentum transfer. Hence, the friction force is stronger provided that the 
wavelength exceeds both of these expressions, i.e., 

A > 



l^b ^ab <^m,ab {Vm,ab Hb ^v,a 



1/2- 



The presented conclusions hold for the initial regime only. The species are coupled through 
collisions and the speed difference Va — vi,, in both the friction force and in the viscosity 
terms containing the speed difference may approach zero provided the time interval is 
long enough. This issue (related to friction force) has been discussed in Vranjes et al. 
[12]. The velocity difference relaxation is altitude dependent, and after presenting the 
detailed collision frequencies in the previous sections, we can make some estimates for 
proton interaction with neutrals. Assuming that protons start to move with some initial 
speed Vpo through static background of hydrogen atoms, we can calculate the time needed 
for both species to achieve some common speed Vc, which naturally must be between 
and VpQ. This is performed for several layers to see the differences caused by density 
and temperature variation. Starting from the momentum equations with friction only 
dvn/dt = Vnpivp — Vn), dvp/dt = Ppnivn — Vp) , wc find the time dependent velocity of the 
two species: 

Vn = Vc-] — — -^ — ^ -exp[-(z/p„ + i/„p)t], (19) 

Vp = Vc — — --^ — ^ -exp[-(z/p„ + z/„p)t]. (20) 

Taking nip ~ rrin and Vno = 0, we obtain for the common velocity for both species 

Vc = = VpO = VpO . (21) 

^pn ~r l^np ^pn ~r l^np ^pO ~r 'T-nO 

In Table fl] the common speed (21) is given for several altitudes, and the time necessary 
for both species to achieve 99 percent of this common speed is calculated from 

ln[0.01] 
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To calculate Upn, which appears here, we used the cross section for momentum transfer 
from Fig. [1} 

From Table [I] we can conclude the following, a) Because of the very frequent collisions 
between the two species, the friction between them should vanish very quickly, but b) at 
the same time this must have a strong effect on the small electrostatic or electromagnetic 
perturbations that are now supposed to set into motion the two species. This implies 
c) that the amplitude of these perturbations must reduce dramatically, d) Because the 
two species move together (after the coUisional time tc), the friction is effectively zero, 
therefore the only relevant remaining dissipation mechanism must be through viscosity, e) 
All these conclusions are altitude dependent; higher up the friction may become dominant 



dissipation process, in particular for a certain wavelength regime, as predicted by Eq. ( 18 ). 
These facts are frequently overlooked in the literature where such a regime of common 
motion is just assumed [but without much regard to the consequence c) above] and equa- 
tions for different components are summed up, and a single fluid dynamics is then studied. 
This implies that the transition process, that takes place within the collision time, is ne- 
glected together with the physical phenomena involved in the process. This problem is 
discussed in Vranjes et al. [12]. A statement of the validity of common equations for all 
species is given also in Alfven and Falthammar [T] on p. 177, where the authors write 
that the common speed has sense only if the speed of neutrals is nearly equal to the 
speed of plasma. Such a situation surely cannot be expected in the initial stadium of 
some electrostatic or electromagnetic perturbations that naturally affect plasma species 
first, with the background of initially immobile neutrals whose dynamics develops due to 
friction and partly due to viscosity, and this after the collisional time only [Vranjes et al. 
[l2]]. In the present work we are able to quantify these phenomena by calculating the 
precise characteristic collisional time and consequently by predicting the amplitude of the 
common speed achieved within such a time interval. 

6.3 Hydrogen dynamics 

We have seen that proton interaction with neutral hydrogen is the most dominant in 
most of the space, so in a proper wave analysis that obeys conservation laws, one needs 
the continuity, momentum, and energy equations describing hydrogen dynamics as well. 
Hydrogen collision frequencies with other species (electrons, ions, and helium atoms) and 
corresponding cross sections can easily be obtained using the momentum conservation 

mjriji'jH = mHnui'Hj- (22) 

For completeness, the collision frequency v^p for the most dominant neutral hydrogen 



atom collisions with protons is presented in Fig. 10 for the elastic scattering (dashed line) 
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and momentum transfer (dotted line). The local minimum in the profile is due to the 
decrease in the profile of target particles (protons). 
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Figure 10: Collision frequency ioi H — p collisions. 

We furthermore need components for hydrogen viscosity and thermal conductivity. 
For rjHH, i^HH we may set 
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With the cross section (Jhh determined by line 2 in Fig. [3} the dynamic viscosity coefficient 
for hydrogen self-collisions becomes 
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(24) 



In writing it we used Uhh = Chh^^h'^t- 

In the literature one can find also the expression that follows from the Chapman and 
Cowling (|9j) model based on the interaction of hard spheres: 
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(25) 



where A/ = l/{\/2'KnHr\) is the mean free path, and r„ is the diameter of the colliding 
particles, in the present case its value is r^ = 2.12 ■ 10^^*^ m. This yields 



QA^mv^ 



(26) 



The expressions (24) and (26) are checked against experimental measurements available 



in Vargaftik et al. [3S] for a hydrogen gas, and the results for several temperatures 
are presented in Table |2} In the given energy range our r^uH gives values closer to the 
experimental ones. The differences between the two models are about factor 2, which 
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Table 2: Hydrogen dynamic viscosity coefficient [in units kg/(sni)= Pa-s] for several 
temperatures. Second row: temperature-dependent cross section from Fig. |3| Third row: 
our value based on Eq. (24). Fourth row: value based on the Chapman and Cowling 



hard sphere model (25). Fifth row: experimental values for pure neutral hydrogen from 
Vargaftik et al. [36j. 



T[K] 



4400 



4990 



6560 



11150 



a^H (a.u.) 45.35 45.3 40 33.5 

3.95 ■ 10^5 4.2 ■ 10-5 5.5 ■ 10"^ 0.86 ■ 10"^ 

10-10-5 11-10-5 13-10-5 1.67-10"^ 

4.8-10-5 5.5-10-5 7-10-5 1.1-10"^ 



Vhh 

Vcc 

Vexp 



is clearly due to indistinguishability effect that is missing in the Chapman and Cowling 
classical model. 

Using the Chapman and Cowling ^ model, we can also calculate the coefficient of 
thermal conductivity for hydrogen gas 
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This can be compared with the above-given corresponding coefficient (23) obtained from 
kinetic theory with the BGK collision integral: 
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Around the temperature minimum region in the photosphere, using line 2 from Fig. [3] for 
hydrogen, our conductivity coefficient k,„h is higher by about a factor 2 than the Chap- 



man and Cowling coefficient (27), which is the consequence of the quantum-mechanical 



indistinguishability incorporated in our derivations. 

For the viscosity coefficients fi^b, which are associated with the terms containing the 
speed difference between H and b species, the situation is as follows: The speed difference 
between different neutral species cannot be of any importance for obvious reasons (they 
are coupled through collisions, they react similarly to perturbations by external forces). 
Hence, both yU^^ and fJ^HHe are not needed for the same reason. We now compare 
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From Fontenla et al. [H] we know that np ^ njj(.+ in the whole region of interest here. 



Therefore fiHHe+ is most likely negligible. Using (|22|) we have 
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Similar arguments are used for the coefficients Xnt in the conductivity vector. Hence, the 
only remaining coefficient we need is 
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Figure 11: Most relevant viscosity {rj„H,^„p) and thermal conductivity {hhh,Xhp) coeffi- 
cients for hydrogen atoms. 



The most important coefficients for hydrogen r]„H, ^Hp-, f^HH, Xnp are presented in Fig. 11 
Evidently, for practical purposes in the lower solar atmosphere the self-interaction coeffi- 
cients TjuH, i^HH may be taken as constant and their values are ~ 0.5 ■ 10"'' kg/(sm) and 
~ 5 ■ lO^'^ (sm)~^, respectively. The other two coefficients /i^p and Xnp, which include 
interaction with other species are strongly altitude dependent. 



7 Summary and discussions 

The parameters in the lower solar atmosphere change with altitude and much care is 
needed to properly describe the physical processes that take place there. One obvious 
example is presented in Fig. [5] where the e-H collision frequency changes by seven orders 
of magnitude between the altitudes of —100 km and 2200 km, taking values 1.3 ■ 10^° 
Hz and 2.4 ■ 10^ Hz, respectively. In addition to this, the type of collisions changes as 
well, e-H collisions being the most dominant up to 900 km and e-p collisions above that. 
A detailed knowledge of these processes is essential to estimate the friction and related 
phenomena (e.g. conductivity, transport, etc.). 

Our most important conclusions can be summarized as follows: 

i) The cross sections presented in Sec. 2.1 are the most accurate existing ones. They 
contain the following essential details: a) variation of cross sections with temper- 
ature (altitude), b) variation of cross sections due to quantum effects in the given 
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low-temperature range in the lower solar atmosphere, and c) clear and pronounced 
differences of cross sections describing elastic scattering, momentum transfer, and 
viscosity. Combined, these fine details may introduce significant differences for 
various processes related to magnetization, transport, heating, etc. Although well- 
known in the laboratory plasmas, so far these details have not been studied in the 
solar atmosphere. Unlike the various approximate data, the data we used here (all 
of them from cited references) are fully quantum-mechanical and obtained with- 
out (almost) any approximation, thus are the most accurate data ever obtained 
in coUisional physics (numerical accuracy to six significant digits, physical accu- 
racy bellow one percent). The only assumption used in their derivation was that 
the electronic excitation to the excited nonresonant electronic states is negligible 
(though charge transfer is included). This approximation is quite accurate below 
the energy threshold for electronic excitation and this is the source of the physical 
accuracy of 'only' one percent. In deriving these data even ro-vibrational degrees 
of freedom were taken, and they are created to serve as a benchmark for checking 
the accuracy of other approximate approaches. Needless to say, we used the exact 
ion-atom potentials from R = to R = 80000 a.u. 

The cross sections presented here, coincide with classical at high energies. At low 
energies (roughly below 1 eV) the cross sections include the effects of a) indistin- 
guishability and b) charge transfer. The lower solar atmosphere is indeed within 
this low-energy range, and consequently these intrinsic properties of the plasma-gas 
matter cannot be avoided. The accuracy of our collision data for ion-atom collisions 
removes possible doubts on the size of the momentum exchange used in previous 
works available in the literature, including possible under- or overestimation of the 
role of Alfen waves and kink wave damping in the lower solar atmosphere. 

ii) For electron dynamics above 850 km neutrals plays no practical role, although the 
neutral number density at 850 km is still three orders of magnitude higher than 
that of protons and electrons. However, this may not be so if inelastic collisions are 
taken into account, e.g., those in which electrons are lost or created, as discussed in 
Sec. 2. These phenomena are beyond the scope of the present work. 

iii) For proton dynamics the role of neutrals is negligible above 1900 km, although at 
this altitude (according to data from Fontenla et al. [H]) Up = 4.24 ■ 10^^ m~^ is 
still below the neutral hydrogen density 12^ = 1.7 ■ 10^"^ m~^. We stress again that 
this conclusion may not hold if inelastic collisions are taken into account. 

iv) There exists a layer within which both electrons and ions are definitely un-magnetized. 
For intense magnetic structures with a magnetic field of 0.1 T this layer is located 
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around an altitude x = and below. In this layer the magnetic field plays no direct 
role in the dynamics of both electrons and ions. 

v) The layer of unmagnetized electrons and ions continues with a much thicker layer in 
which electrons are magnetized and ions are not. The depth of this layer changes 
spatially, and for protons it is at least 1000 km thick in regions with a kilo-Gauss 
magnetic field. The dynamics of electrons and ions in this region is completely dif- 
ferent, and models that assume a two-component system consisting of 'neutrals' on 
one side and 'plasma' on the other are meaningless. This is because the 'plasma' 
contains electrons and ions whose dynamics is totally different because of the mag- 
netic field. Consequently, in this layer the electrons and ions cannot be treated as 
a single fiuid. A fully multi-component analysis (fiuid or kinetic) has no alternative 
in the lower solar atmosphere. 

vi) Viscosity and thermal conductivity coefficients given in this work are currently the 
most accurate, and at the same time the most complete ones because they con- 
tain all most relevant terms appropriate for a multi-component system such as is 
the lower solar atmosphere with un-magnetized ions. They also completely agree 
with experimental measurements. Our results show that including viscosity both 
for protons and neutral hydrogen may be essential to properly capture diffusion 
phenomena in the solar atmosphere. 

If we return back to Fig. [T| we see that the proton-hydrogen cross section changes with 
temperature (i.e. with the altitude). Its values for elastic scattering (in the laboratory 
frame) are 2.269 ■ 10"^^ m^ and 1.601 ■ 10"^*^ m^ at temperatures of 5 ■ 10^ K and 20 • 10^ K, 
respectively. The temperature of 5 ■ 10^ K corresponds to the altitude 200 km (and also 
705 km), and 20 ■ 10^ K corresponds to the altitude ~ 2200 km. We can can now compare 
these with the cross section used by other researchers, e.g., Zaqarashvili et al. [H], where 
the cross section was assumed to be constant with the value 8.79 ■ 10~^^ m^, i.e., vr a.u., 
describing collisions where ions and neutrals are treated as hard spheres. Our correct 
values, which include the quantum-mechanical effect of indistinguishability, for the two 
temperatures given above are 258 and 182 times greater! We note that these authors 
make no distinction between cross sections for elastic scattering and momentum transfer. 
Therefore we can also compare their value with our cross section for the momentum 
transfer from the same figure; for the two temperatures we have 1.040 ■ 10~^^ m^ and 
8.785 • 10^^^ m^^. These are again 118 and 100 times greater than their values. We 
observe that in the mentioned work the proton-helium cross section is assumed to be the 
same as the one for proton-hydrogen collisions given above. However, from our Fig. [2] 
the cross sections for p-He elastic scattering at 5 ■ 10^ K and 20 ■ 10^ K in the laboratory 
frame are about 110 and 70 times greater than their value. At the same time, our cross 
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section for momentum transfer for the two temperatures is 20 and 8 times greater than 
their value. 

On the other hand, the cross section for p-H colhsions in Khomenko and Collados [2T] 
is fixed to 5 ■ 10~^^ m^, which is rather close to our value for the momentum transfer 
cross section in plasma reference frame, roughly speaking only twice as smaller. Though 
in their subsequent calculation of the collision frequency this difference is compensated 
by the numerical factor 8^'^, which they keep in the thermal velocity, and the collision 
frequency which they obtain is very close to our value. 

In view of the results presented here a natural next step is to include effects of inelastic 
collisions. In our previous work (Vranjes and Poedts |39]) we have shown that in certain 
layers in the photosphere all ions in a unit volume recombine at least 26 times per second. 
This may have consequences on magnetization, for example. Perhaps this may be used 
also to explain the nature and longevity of prominences, which are believed to contain 
considerable amounts of neutrals. Their longevity is a challenge for the theory because 
neutrals should naturally diffuse and evacuate a prominence by moving to lower layers 
due to gravity. However, in the presence of inelastic collisions this diffusion should take 
place at reduced speed because a neutral particle does not remain neutral all the time, it 
is consequently affected by the magnetic field and the prominence may last longer. 
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